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Section 1 
INTRODUCTION 


At the request of the Governments of Canada and the United States, the 
International Joint Commission has undertaken an investigation to determine 
the possibilities for limited regulation of Lake Erie and the consequent 
effects throughout the Great Lakes basin. As a result of this request, the 
Commission established the International Lake Erie Regulation Study Board on 3 
May 1977 to undertake, through appropriate agencies in Canada and the United 
States, the necessary investigations and studies and to advise the Commission 
on all matters which it must consider in reporting on this matter. 


The Working Committee, which was established by the Study Board, appointed 
a Regulation Subcommittee composed of members for the U.S. from the Buffalo 
and the Detroit Districts of the Corps of Engineers and members for Canada 
from the Inland Waters Directorate, Environment Canada, located in Cornwall 
and Burlington, Ontario. The Regulation Subcommittee and, in particular, 
these two agencies were assigned the task of deriving the basis-of-comparison 
and the data required for the development of regulation plans. Water supply 
data and the basis-of-comparison data, covering the period January 1900 
through December 1976, have been developed and coordinated. This report 
provides a description of the methods employed to obtain the final data, and 
contains tabulations of the basis-of-comparison data as well as the basic data 
used in their derivation. 





Section 2 
RECORDED DATA 


The historic or recorded data, from which the required water supplies and 
comparison data for the regulation plans have been developed, were themselves 
developed, in part, by another international committee known as the 
Coordinating Committee on Great Lakes Basic Hydraulic and Hydrologic Data. 
This Committee is comprised of three Canadian and three U.S. members who 
represent the federal agencies in each country having responsibilities for 
accurate determination of such data. The Committee develops and coordinates 
data on water levels, river flows and on the physical characteristics of the 
Great Lakes system. The values developed by the Coordinating Committee on 
Great Lakes Basic Hydraulic and Hydrologic Data were employed where possible. 
Where required coordinated data did not exist, the missing data required were 
developed by the Board’s Regulation Subcommittee. 


Due to their larger lake area, the levels of Lake Superior and Lakes 
Michigan-Huron respond to changes in water supply much more slowly than do the 
levels of Lakes Erie and Ontario (Lake Ontario six times as fast as Lakes 
Michigan-Huron). For this reason, the basic data used in this study were 
developed and coordinated for monthly periods on Lakes Superior and 
Michigan-Huron, quarter-monthly periods on Lakes Erie and Ontario. Since Lake 
St. Clair reflects conditions mainly on Lakes Michigan-Huron, monthly periods 
were employed on that lake. In order to measure the effect of upstream 
regulation on the lower St. Lawrence River, data were also developed for Lake 
St. Louis. 


The quarter-monthly periods used consist of seven or eight days 
subdivided as shown in the following tabulation: 


Dates Included in Quarter-Monthly Periods 


Months of Months of 
Quarter Month 28 or 29 days 30_or 31 days 


First 1-7 1-8 
Second 8-14 9-15 
Third 15-21 16-23 
Fourth 22-28 or 29 24-30 or 31 


The recorded data, which are contained tn Annex A and which were employed 
in the computation of water supplies and the basis-of-comparison for the 
period January 1900 through December 1976, are discussed below. 

2.1 Beginning of Period Lake Levels 


Since the beginning-of-period lake level is theoretically an 
tnstantaneous value and because of daily and even hourly level fluctuations 





resulting from unbalance or tilting of the lake surface caused by winds, 
differential barometric pressure and seiche effects, the level selected might 
be subject to considerable error if only one gauge on a large Jake were used. 
Therefore, in determining the beginning-of-period levels, several gauges were 
selected at points around the lake to provide a good representation of the 
lake level. The beginning-of-period levels were selected and coordinated by 
graphically, mathematically or judgementally eliminating the effects of wind, 
datum differences and other factors. 


Basically, beginning-of-period levels were computed using the following 
method. When only one gauge was used to determine the beginning-of-period 
level for a lake, the level was calculated as the average of the daily mean 
levels for the first two days of the period and the last two days of the 
previous period. Where more than one gauge was used, the level was computed 
by first averaging the daily mean levels at all the gauges on the first day of 
the period, and then on the last day of the previous period. These two values 
were then in turn meaned to obtain a beginning-of-period level for the lake. 


For Lake Superior, the water level gauges at Thunder Bay, Duluth, 
Michipicoten, Marquette, and Pte. Iroquois were used. The gauges used on 
Lakes Michigan-Huron were located at Milwaukee, Ludington, Mackinaw City, 
Harbor Beach, Thessalon, and Goderich. 


For Lake St. Clair the Windmil? Point gauge was used for the period 
1900-1911, the Isle Aux Peche gauge from 1912 to 1925, the Tecumseh gauge from 
1926-1950, the Grosse Pte. gauge for the period 1951-1969, and the St. Clair 
Shores gauge for the period 1970-1976. Relationships between each of these 
gauges and the Fort Wayne and Port Lambton river gauges were used to refine 
the levels when recorded levels did not appear representative of river and 
lake level conditions. 


Lake Erie beginning-of-period levels were determined from Cleveland and 
Port Stanley daily records. The gauges at Toledo and Buffalo were used as 
supporting data during unstable lake level conditions. Lake Ontario 
beginning-of-period levels are available on a current basis from the 
International St. Lawrence River Board of Control. 


2.2 Recorded Out flows 


Recorded outflows from Lake Superior via the St. Marys River, the St. 
Clair and Detroit River flows and the Lake Erie and Ontario outflows were 
supplied by the Coordinating Committee. Lake St. Louis outflows were compiled 
by the Inland Waters Branch and adopted by the Regulation Subcommittee. 
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2.3 Recorded Diversions 


The Long Lac diversion was begun in 1939. Water from Long Lac, which 
naturally drained into James Bay via the Kenogami and Albany Rivers, is 
diverted through a series of small lakes into the Aguasabon River, a tributary 
to Lake Superior. The Ogoki diversion was begun in 1943. Water from the 
Ogoki River, another tributary of the Albany River, is diverted into Lake 
Nipigon and held there until required by hydro-electric plants on the Nipigon 
River which drains into Lake Superior. Data on the amount of water diverted 
at both Long Lac and Ogoki are supplied by the Hydro-Electric Power Commission 
of Ontario and are reported in the International Lake Superior Board of 
Control Annual Report. Data from this report was used where possible. 


The Chicago Sanitary and Ship Canal, between the Chicago River and the 
Des Plaines River, forms part of the Illinois Waterway connecting Lake 
Michigan and the Mississippi River. The flow in this canal is controlled by a 
dam and gates at Lockport, Illinois. Discharge figures are supplied to the 
U.S. Army, Corps of Engineers by the Chicago Sanitary District and have been 
adopted for use in this study. 


The Welland Ship Canal crosses the Niagara Peninsula from Port Colborne 
on Lake Erie to Port Weller on Lake Ontario. The New York State Barge Canal 
System is a series of canals, a part of which links the Niagara River above 
Niagara Falls with Lake Ontario. Recorded discharges for the Welland Canal 
and New York State Barge Canal for the period 1900-1975, are available in the 
Coordinating Committee's Lake Erie Outflow Report dated June 1976. This 
report was used where possible. The Monthly Report of Flow Directed from Lake 
Erie-Welland Canal, released by the St. Lawrence Seaway Authority of Transport 
Canada, was the source of the 1976 values. 


The recorded data is tabulated in Annex A of this report. 
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Section 3 
DERIVED DATA 


Employing the recorded data listed above, data required for the lake 
regulation studies were derived as follows and are tabulated in Annex B. 


3.1 Net Basin Supplies 


The net basin supply is a term used to describe the water which a lake 
receives from precipitation on both its surface and on its own land drainage 
basin less the net effect of evaporation and condensation on the lake surface. 
With presently available techniques, some of these factors cannot be 
determined accurately; therefore, the net basin supplies were computed by 
employing reliable lake level and flow records for the required monthly and 
quarter-monthly periods. The relationship used is as follows: 


NBS = AS+#0-140D 


Where: 
NBS = Net basin supply 
AS =.Change in storage from beginning to end of the period 
0 = Average outflow from lake through outflow river 
I = Average inflow to lake from inflow river 
D = Diversion into (-) and out of (+) the lake 


with all terms expressed in units of cubic feet per second for the period. 
The sum of the values AS and 0 represents the net total supply to the lake. 


The values for changes in storage, outflows and inflows are obtained or 
determined directly from the mutually adopted recorded values of lake levels, 
river flows and diversions listed above. 


3.2 Winter and Weed Retardation 


The flows in the outlet rivers of the lakes during the winter season are 
often retarded materially by ice formation and by ice jamming. These 
conditions are not predictable for any specific winter, either as to their 
severity or the exact timing of their occurrence. The natural retardation of 
flows under ice conditions causes the levels of unregulation lakes to be 
higher at the time of the spring breakup than the levels would be if there is 
no ice, and this increases the storage on the lake. 
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The water level of Lake Superior and the outflow through the St. Marys 
River are regulated by the International Lake Superior Board of Control under 
authority of the International Joint Commission. The physical control is 
achieved by a dam and other structures at the head of the St. Marys Rapids at 
Sault Ste. Marie, Michigan/Ontario. Under present regulation conditions, the 
winter retardation effect on the discharges is virtually zero. Since the 
basis-of-comparison condition for Lake Superior was considered to be under 
regulated conditions, it was not necessary to consider winter retardation in 
the St. Marys River. 


Lake Huron does not freeze over completely during the winter, primarily 
due to the influence of wind and the heat stored in the lake. The ice which 
forms in exposed central parts of the lake is continually broken up and moved 
about by the action of the wind. Some of this ice finds its way into the St. 
Clair River. As a result of these heavy runs of ice, jams occur which 
materially reduce the normal flow and in turn affect both the upstream and 
downstream levels. Winter retardation in the St. Clair River was computed for 
use in this study by subtracting the coordinated recorded St. Clair River flow 
from the corresponding discharge computed from the open water discharge 
relationship for the Harbor Beach, and St. Clair Shores (or Grosse Pointe) 
gauges. 


Lake St. Clair normally freezes over in early winter and shields the 
Detroit River from heavy ice runs. The Detroit River itself, however, 
frequently freezes over in its lower reaches. Due to the size of Lake St. 
Clair, even a small retardation influences its level regime. Therefore, for 
use in this study, winter retardation in the Detroit River was determined to 
be the difference between the flow computed from the open water discharge 
relationship for the St. Clair Shores (or Grosse Pointe) and Cleveland gauges, 
and the coordinated recorded flow. 


As in the case of Lake Huron, the principle problem with ice in Lake Erie 
and the Niagara River results from breakup of lake ice fields and the wind 
pushing the ice into the river. In order to avoid ice problems in the Niagara 
River, an ice boom has been placed across the head of the Niagara River each 
winter season since 1964 by the power entities. The presence of this ice boom 
has reduced the retardation of Niagara River flow by ice to a very small 
amount which can be considered insignificant. However, since the outlet 
conditions of 1953 were adopted as the basis to be used for comparing 
regulation plans, average winter retardation was assumed for the Niagara River 
over the period of record. In addition, a small amount of weed retardation 
exists in the summer months. For these investigations, average monthly values 
as derived by the Coordinating Committee on Great Lakes Basic Hydraulic and 
Hydrologic Data in the Lake Erie Outflow Report dated June 1965 have been 
used. 


Since Lake Ontario has been regulated since 1960, and the 
basis-of-comparison assumes this condition for the entire study period, no 
winter retardations were required for the calculations of effects on the upper 
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St. Lawrence River. Reduction in the winter flow at the outlet of Lake St. 
Louis was calculated directly as the difference between the discharges derived 
from its approximate open-water stage discharge curve and the recorded 
discharges. 
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Section 4 
BASIS FOR COMPARISON 


The recorded Great Lakes level and outflow data reflect the effects of 
changes in the regimen of the lakes and their connecting channels which 
occurred over the study period (1900-1976). The principal changes were 
man-made and consisted of changes in the amount of diversion into and out of 
the Great Lakes basin, alterations in the configuration of the connecting 
channels and the erection of control structures at the outlets of Lakes 
Superior and Ontario. In addition to these man-made changes, the movement of 
the earth's crust in the Great Lakes basin, a natural phenomenon, has been 
progressively changing the tilt of the basin with a resultant gradual change 
in the elevation of individual lake outlets with respect to any given point on 
each lake. In order to permit consistent evaluation of the effects of various 
regulation plans, it was necessary to develop a set of Jake levels and 
outflows which would reflect a constant or fixed regimen in the Great 
Lakes-St. Lawrence River System over the study period. 


The Board agreed on the following basic conditions of lake levels and 
outflows to be used as the basis for evaluation. The conditions, together 
with the reasons for their selection are as follows: 


(1) A constant diversion of 5,000 cfs into Lake Superior by way of the 
Long Lac and Ogoki diversions. This diversion was authorized under the 
exchange of notes, dated October 14 and 31 and November 7, 1940, between the 
United States and Canada and has averaged approximately this amount since that 
date. 


(2) Lake Superior regulated in accordance with Plan 1977, which is the 
currently authorized plan being used by the International Lake Superior Board 
of Control for determining releases from Lake Superior. 


(3) A constant diversion of 3,200 cfs out of Lake Michigan at Chicago. 
This is the maximum allowable diversion at Chicago by decree of the U.S. 
Supreme Court, dated June 12, 1967. 


(4) 1962 outlet conditions for Lake Huron. This represents the current 
conditions, which have existed since the completion of the 27-foot navigation 
channel dredging in 1962. 


(5) A constant diversion, by way of the Welland Canal, of 7,000 cfs out 
of Lake Erie and into Lake Ontario. This is the approximate average diversion 
through the Welland Canal during the period 1950-1976. 


(6) 1953 outlet conditions for Lake Erie. In its 1953 report on the 
Preservation and Enhancement of Niagara Falls, the International Joint 
Commission considered it essential that the relationship existing at that time 
between the Niagara River flow and the Chippewa-Grass Island Pool level be 


A-8 


Apt amt settee s Men See eer” - * - CO a TPT eS “s - -- wt ee ee ee 


weiner ee 


ce ae 7 sain a ae a Saath ln, 


mere, emt er cee RE npn ng 


we - 


maintained following the commencement of operation of the Chippewa-Grass 
Istand Pool Control Structure and power diversions as permitted by the 1950 
Niagara Treaty. 


(7) Lake Ontario regulated during the period 1900-April 1960 in 
accordance with Plan 1958-D without discretionary deviation. For the period 
from April 1960 to December 1976, Lake Ontario was regulated in accordance 
with Plan 1958-D with discretionary deviations as they occurred. Minor 
adjustments to the discretionary deviation values were required during the 
period 1973-1976 to preclude violation of the St. Lawrence River low water 
profiles. 


(8) Recorded conditions for the Ottawa River and local inflow to the St. 
Lawrence River. 


The levels and outflows to be used as a basis-of-comparison for each lake 
were obtained by routing through the system the coordinated net basin supplies 
employing the constant conditions previously listed. 


For Lake Superior, the basis-of-comparison levels and outflows were 
obtained by routing through the lake the coordinated net basin supplies 
(adjusted for a constant 5,000 cfs diversion into the lake by way of the Long 
Lac and Ogoki diversions) in accordance with the present regulation plan known 
as Plan 1977. 


Because of the nature of the control in the St. Clair and Detroit Rivers, 
the problem of routing supplies through Lakes Michigan-Huron, St. Clair and 
Erie was fairly complex. The St. Clair River and Detroit River flows are 
dependent not only on Lake Huron levels, but also on the water levels in the 
lower river and Lake Erie. Therefore, a method of successive approximations 
of supply routing as described in the steps below was used: 


(1) The month was first considered to be broken-up into very small 
periods (the period used was one-fortieth of a month) because Lake St. Clair 
levels are very sensitive to small changes in outflow. Therefore, the data 
available on a monthly basis (net basin supplies to Lakes Michigan-Huron and 
St. Clair) were considered to be made up of 40 equal parts, and 
correspondingly, the data available on a quarter-monthly basis (net basin 
supplies to Lake Erie) were considered to be made up of 10 equal parts. 


(2) Net total supply to Lakes Michigan-Huron was calculated as the 
routed Lake Superior outflow plus the Lakes Michigan-Huron net basin supply 
minus a constant diversion of 3,200 cfs at Chicago. 


(3) In the first approximation, the beginning-of-period stages were 
considered equal to the average levels for the period. 


(4) A first estimate of the Lake Huron outflow was made with the Harbor 
Beach-St. Clair Shores two-gauge relationship for the 1962 regime: 
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Q = 0.0841168 (0.5 HB + 0.5 SCS - 543.4)2-0 (HB - scs)0-5 - winter 
retardation. ; 











(5) A first estimate of the Lake St. Clair total supply was made by 
adding the Lakes Michigan-Huron calculated outflow to the Lake St. Clair net 
basin supply. 


(6) A first estimate of the Lake St. Clair outflow was made from the St. 
Clair Shores-Cleveland two-gauge relationship for the 1962 regime: 











Q = 0.1280849 (SCS - 543.81)2 (SCS - C)9°5 - winter retardation. 





(7) An initial estimate of the Lake Erie total supply was made by adding 
the Lake St. Clair calculated outflow to the Lake Erie net basin supply. 







(8) A first estimate of Lake Erie outflow was made from the Buffalo 
single gauge relationship for the 1953 regime. 






Q = 3665 (B - 556.25)1-5 + 7.00 (Welland) - winter or weed retardation. 






(9) Having a first estimate of the total supply and outflow from each 
lake, a first estimate was made of the change of storage, the end-of-pertod 
level, and the mean level for the period. 







(10) Steps (4-9) were repeated using first estimate mean levels in place 
of beginning-of-period levels. 





(11) The process was repeated four times. This was judged sufficient to 
stabilize the system. 






(12) The required monthly and quarter-monthly data were extracted by 
combining the calculated data for the shorter periods. 












Having the Lake Erie outflows and Lake Ontario net basin supplies, the 
net total supplies to Lake Ontario were calculated. OQutflows and levels for 
Lake Ontario were calculated by employing the developed rules of Plan 1958-0 
for the period up to 20 April 1960 and applying the deviations from this plan 
as they actually occurred during the period after that date to December 1976, 
with some minor adjustments made to preclude violation of the St. Lawrence 
River low water profiles during the latter period. 















Lake Ontario outflows were added to the recorded differences between the 
Lake St. Louis and Lake Ontario outflows to obtain Lake St. Louis outflows for 





comparison purposes. Then, working with the Lock 5 rating equation, the 
basis-of-comparison Lake St. Louis levels were calculated as: 


SL = [(Q + ice retardation) /5.6268]¢/3 + 54.29 


The basis-of-comparison data are contained in Annex C. 
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Section 5 
SUMMARY 


The current study of possible regulation of Lake Erie has required the 
coordination of the net basin supplies and basis-of-comparison levels and 
outflows for each of the lakes. These basic data are employed directly in the 
regulation study and serve a number of important functions in the study. 


Most of the long term trends, for example, changes in diversions and 
regulation procedures, have been removed from the data. Several other long 
term trends such as increasing tributary regulation, changes in winter 
retardation, and increasing consumptive use remain in the data. However, it 
is assumed that these effects are small enough so as not to significantly 
affect the regulation studies over the historical period, when compared with a 
basis-of-comparison which employed the same input data. 
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ANNEX A 
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ANNEX A 
RECORDED DATA 


Lake Superior Monthly Beginning-of-Period Levels 

Lakes Michigan-Huron Monthly Beginning-of-Period Levels 
Lake St. Clair Monthly Beginning-of-Period Levels 

Lake Erie Quarter-Monthly Beginning-of-Period Levels 
Lake Ontario Quarter-Monthly Beginning-of-Period Levels 
Lake Superior Monthly Mean Outflow 

Lakes Michigan-Huron Monthly Mean Outflow 

Lake St. Clair Monthly Mean Outflow 

Lake Erie Quarter-Monthly Mean Outflow 

Lake Erie Monthly Mean Outflow 

Lake Ontario Quarter-Monthly Mean Outflow 

Lake Ontario Monthly Mean Outflow 

Lake St. Louis Monthly Mean Outflow 

Long Lac and Ogoki Diversions 


Chicago Sanitary and Ship Canal (Including I}]linois and 
Michigan Canal) 


Welland Canal 
New York State Barge Canal 
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Section 1 
ABSTRACT 


The Great Lakes Regulation Model was developed by the Great Lakes 
Hydraulics and Hydrology Branch, Lake Hydrology Section, for use on a C.D.C. 
7600, 6600 or 6400 computer system. 


This model has been designed for multiple-uses, incorporating three basic 
programs to represent the manner in which water transcends through the Great 
Lakes system, from Lake Superior to the St. Lawrence River. Two of the 
programs mode? the existing regulatory controls on Lakes Superior and Ontario 
by Regulation Plans 1977 and 1958-D, respectively. A third program routes the 
supplies and outflows through the middle lakes using the 1962 - 1968 
conditions on the St. Clair and Detroit Rivers and the 1935 conditions on the 
Niagara River. 


Each of the basic programs give the user the capacity to vary inputs for 
studies of the Great Lakes. In particular, the middle lake routing program 
allows for varying the Lake Erie outflows by using the net basin supplies for 
Lake Erie and/or Michigan-Huron and/or Superior. It is this portion of the 
overal] model which is utilized in the International Lake Erie Regulation 
Study. 


The Great Lakes Regulation Model is designated by the Detroit District as 
program number 724F3058. It requires the following subroutines be run: 1 
QSUP; 2) ABSICE; 3) CANADP; 4) CANADA; 5) R1958D; 6) STDMY; 7) PTITLE; 8) 
IABICE; 9) PRTHCE; 10) PRTONT; and 11) ROUND. 


The printed output contains the levels and flows data for the Great Lakes. 
The tape output consists of 19 files which contain the levels and flows of the 
Great Lakes, formatted to include the beginning-of-month and monthly mean lake 
levels and the monthly mean outflows from the Great Lakes. For the lower 
lakes (Lakes Erie and Ontario) this information is presented in quarter-month 
format. There are also output files containing monthly mean water levels for 
the U.S. Slip Gauge on the lower St. Marys River, the projected gate settings 
for the Lake Superior Compensating Works, and the quarter-monthly mean water 
levels and outflows of Lake St. Louis. 








Section 2 
SPECIFICATIONS 


2.1 Limitations 

The number of years of routing are limited only by the number of years of 
water supply data inputted. Subroutines "ROUND" and “STDMY" depend upon the 
number of significant digits of the machine and the system clock routines, 
respectively. 
2.2 Input Data 

The program reads the basic data from Tape File 66. The controlling 
parameters, along with the initial data, are read from the input file. This 
data has been tabulated in Volume 2 of this Appendix. 

2.2.1 Tape File 66 


This file consists of 25 card images per year. The format is (9X,14,7X, 





1215/(20X,1215)) with the second part of the format used for reading 


quarter-month data. 


Each year has the following forn. 


Card Column Column Content 
10-13 21-80 

1 IYR-(Year) NBS (J,1), NBSS Lake Superior net basin 
supplies (HCFS). 

2 " ICE (J,1) Lake Superior ice 
retardation (HCFS). 

3 " NBS (J,2), NBSM Lakes Michigan-Huron net 
basin supplies (HCFS). 

4 , ICE (J,2), ICEM Lakes Michigan-Huron ice 
retardation (HCFS). 

5 " NBS (J,3), NBSSC Lake St. Clair net basin 
supplies (HCFS). 

6 . ICE (0,3), ICEC Lake St. Clair ice 
retardation (HCFS). 

7 " NBSEQ, NBSE Lake Erie net basin 
supplies-Jan, Feb, Mar 
(HCFS). 

8 ° NBSEQ, NBSE Lake Erie net basin 
supplies-Apr, May, 
June (HCFS). 

9 " NBSEQ, NBSE Lake Erie net basin 
supplies-July, Aug, 
Sep (HCFS). 

10 * NBSEQ, NBSE Lake Erie net basin 
supplies-Oct, Nov, Dec 
(HCFS). 

11 * ICE (J,4) Lake Erie flow retardation 


(HCFS). 








Read in Subroutine “R1958-D" 









Column 
10-13 


12 IYR-(Year) NBS Lake Ontario net basin 
supplies-Jan, Feb, Mar 
(HCFS). 
13 . NBS Lake Ontario net basin 
supplies-Apr, May, June 1 
(HCFS). 
14 ' NBS Lake Ontario net basin 
supplies-July, Aug, Sep 
(HCFS). 
15 . NBS Lake Ontario net basin 
supplies-Oct, Nov, Dec 
(HCFS). 
16 . ICE Lake Ontario ice 
i retardation (HCFS). 
17 . IFLO Lake St. Louis - Lake 
Ontario flow-Jan, Feb, 
Mar (TCFS). 
18 " IFLO Lake St. Louis - Lake 
Ontario flow-Apr, May, 
June (TCFS). 
f 19 * IFLO Lake St. Louis - Lake 
Ontario flow-July, Aug, 
Sep (TCFS). 
20 ? IFLO Lake St. Louis - Lake 
Ontario flow-Oct, Nov, 
Dec (TCFS). 





Column 21-80 Content 










ee me ne ee ene ee ee eae 














21 Z IDIS Deviation from Rule 1958-D 
Jan, Feb, Mar (TCFS). 
22 . IDIS Deviation from Rute 1958-D 
Apr, May, June (TCFS). 
23 . IDIS Deviation from Rule 1958-D 
24 . IDIS Deviation from Rule 1958-D 
Oct, Nov, Dec (TCFS). 
25 . ICESL Lake St. Louis ice 






retardation (HCFS). 







2.2.2 Input File 






The program reads from 1-10 title cards, in the format (26A3,A2), into 
PTITLE. These title cards must be followed by a blank card, unless there are 
The rest of the data cards are as follows. 









10 cards. 





July, Aug, Sep (TCFS). 


Card 


1 


Columns 


1-72 
1-72 
1-72 
1-72 
1-72 
1-72 
1-5 


6-10 
11-15 


16-20 


Format 
12F6.0 
12F6.0 
12F6.0 
12F6.0 
12F6.0 
12F6.0 
15 


15 
15 


15 


Variable 


Name 


PSM 
PSSD 
PHM 
PHSD 
QB 
SUPL IM 
IBY 


IEY 
IFO 


ICES 


SEOP 
CHNG 
AVAR 
TEMPSO 


EICE 


Variable Description 


Lake Superior target stages 
(feet-IGLD). 

Lake Superior standard deviations 
(feet). 

Lakes Michigan-Huron target stages 
(feet-IGLD). 

Lakes Michigan-Huron standard 
deviations (feet). 

Lake Superior target flows (TCFS). 

Lake Superior lower flow limit (TCFS). 

Beginning year. 

Ending year. 

Deviation code for Rule of 1958-D: 
l-for deviations in basic data; 
2-for no deviations. 

l-for ice retardation in St. Marys 
River. 

0-for no ice retardation in St. Marys 
River. 

l-for ice retardation in St. Clair 
River. 

0-for no ice retardation in St. Clair 
River. 

l-for ice retardation in Detroit 
River. 

0-for no ice retardation in Detroit 
River. 

l-for flow retardation in Niagara 
River. 

O-for no flow retardation in Niagara 
River. 

10000-for single winter gate setting. 

IMMNN-for single winter gate setting 
and single gate setting from month 
MM to month NN. 

OMMNN-for a single gate setting fron 
month MM to month NN. 

00000-for the program to change gate 
setting as needed. 

Lake Superior beginning-of-period 
stage (feet-IGLD). 

Al lowable change in Lake Superior 
outflow per month (TCFS). 

Constant (usually 200 TCFS). 

Lake Superior monthly mean flow for 
the month previous to the starting 
month (TCFS). 

Niagara River ice and weed flaw 
retardations (HCFS). 





Card 


10 


Columns 


1-5 
6-10 


Variable 


Name 


LLO 
LNM 


Variable Description 


Long Lac-Ogoki Diversion reduction 
amount (HCFS). 

l-for Lakes Michigan-Huron net basin 
supplies being used in NBH. 

O-for Lakes Michigan-Huron net basin 
supplies not being used in NBH. 
l-for Lake Superior net basin supplies 

being used in NBH. 

O-for Lake Superior net basin supplies 
not being used in NBH. 

Running 12-month sum of Lake Superior 
plus Lakes Michigan-Huron net basin 
supplies, which must equal or 
exceed to flow less water fram Long 
Lac-Ogoki Diversion. 

Lakes Michigan-Huron 
beginning-of-month water level 
(feet-IGD). 

Lake St. Clair beginning-of-month 
water level (feet-IG@D). 

Lake Erie beginning-of-month water 
level (feet-IG@D). 

Basic Chicago Diversion outflow for 
plan (TCFS). 

Net basin supplies for Lake Superior 
for year previous to beginning year 
(HCFS). 

Net basin supplies for Lakes 
Michigan-Huron for year previous to 
beginning year (HCFS). 

Net basin supplies for Lake Erie for 
year previous to beginning year 
(HCFS). 

l-for Lakes Michigan-Huron net basin 
supplies being used in NBH2. 

0-for Lakes Michigan-Huron net basin 
Supplies not being used in NBH2. 

1-for Lake Superfor net basin supplies 
being used in NBH2. 

0-for Lake Superior net basin supplies 
not being used in NBH2. 

Numbur which running 12-month sum of 
Lake Superior and/or Lakes 
Michigan-Huron net basin supplies 
must equal or exceed to flow 
additional water through Chicago 
Diversion. 


















































Card 


15 
Cont. 


16 
17 


18 


19 


20 


Columns 


16-20 


21-25 


1-72 
1-72 


1-6 


1-72 


1-5 


11-15 


16-20 


Format 


I5 


12F6.0 
12F6.0 


12F6.0 


15 


15 


I5 


15 


Variable 


Name 


NBE3 


Kl 


CD 
TBOMS 


PMMS 


NBH 


NBE2 


INA 


INS 


Variable Description 


Number which running 12-month sum of 
Lake Erie net basin supplies must 
equal or exceed to flow additional 
water through the Chicago 
Diversion. 

Number of the first month of the 
“accounting year" for Chicago 
Diversion, l=Jan, 2=Feb, etc. 

(i.e. first month of 12 months used 
for 12-month average 
beginning-of-month stage of 

Lakes Michigan-Huron used for 
flowing additional water through 
the Chicago Diversion). 

Additional Chicago Diversion flow by 
month (TCFS). 

Monthly values that beginning-of-month 
level for Lakes Michigan-Huron must 
exceed for additional Chicago 
Diversion outflow (feet). 

Amount which the average Lakes 
Michigan-Huron monthly mean level 
for the previous “accounting year" 
must exceed or equal to flow 
additional water through the 
Chicago Diversion in the current 
“accounting year" (feet-IGLD). 

Monthly mean Lakes Michigan-Huron 
levels for the year previous to the 
starting year (feet-IG@D). 

Amount which the running 12-month sum 
of Lake Superior and/or Lakes 
Michigan-Huron net basin supplies 
must equal or exceed to flaw 
additional water from Lake Erie 
(HCFS). 

Amount which the running 12-month sum 
of Lake Erie net basin supplies 
must equal or exceed in order to 
flow additional water fran Lake 
Erie (HCFS). 

1-for Lakes Michigan-Huron net basin 
supplies being used in NBH. 

0-for Lakes Michigan-Huron net basin 
supplies not being used in NBH. 

1-for Lake Superior net basin supplies 
being used in NBH. 

0-for Lake Superior net basin supplies 

not being used in NBH. 


we renee meee. 


Card Columns 


20 21-25 
Cont. 


26-30 


31-35 


36-40 


21 1-72 
22 1-72 


23-25 1-80 
26-28 1-80 
29 1-80 
30 1-66 


31 1-5 


6-10 
11-15 
16-20 
21-25 
26-30 
31-35 
36-40 
41-45 


46-50 
51-55 
56-60 


61-63 





Format 


15 


15 


15 


15 


12F6.0 
12F6.0 
1615 
1615 


1615 
22F3 


I3 


NBHE 


NBEE 


INME 


INSE 


ADD 
SUB 
MAX 
MIN 


MLIM 
ISI 


1Q 
ISTAGE 


INDEX 





Variable 
Name 





Variable Description 


Amount which running 12-month sum of 
Lake Superior and/or Lakes 
Michigan-Huron net basin supplies 
must be less than in order to flow 
less water from Lake Erie (HCFS). 

Amount which running 12-month sum of 
Lake Erie net basin supplies must 
be less than in order to flow less 
water from Lake Erie. 

l-for Lakes Michigan-Huron net basin 
supplies being used in NBHE. 

0-for Lakes Michigan-Huron net basin 
supplies not being used in NBHE. 

l-for Lake Superior net basin supplies 
being used in NBHE. 

0-for Lake Superior net basin supplies 
not being used in NBHE. 

Amount by which Lake Erie outflow is 
to be increased (TCFS). 

Amount by which Lake Erie outflow is 
to be decreased (TCFS). 

48 values of maximum "P". 

48 values of minimum “P". 

12 values of minimum “M". 

Supply indicator for 22 quarter-month 
preceeding initial quarter-month 
(TCFS). 

Initial quarter-month of starting 
month (i.e. 1,2,3,or 4). 

Initial month. 

Initial year. 

Final quarter-month of final month. 

Final month. 

Final year. 

Initial weighed supply. 

Exceedence probability. 

Initial quarter-month net total supply 
(TCFS). 

Initial quarter-month supply indicator 
adjustment. 

Initial quarter-month Basic Rule Curve 


flo. 

Initial beginning-of-quarter-month 
stage (hundreds of feet with the 
initial 240 feet deleted, 

i.e., 245.67 is 00567). 
l-for I Limitation. 






























































































































Card 


31 
Cont. 


32 


33 


Columns 


64-66 
67-69 


70-72 


73-75 
76-80 
1-5 


11-15 


16-20 


21-25 


26-30 


1-60 
61-65 


Format 


I3 
13 


13 
13 
I5 


15 


I5 


15 


15 


I5 


1215 
15 


Variable 


Name 


LIM 
IE 


LIM1 


IEl 
TACUM 
IML 


IMH 


TEL 


IEH 


TOL 


IOH 


IUP 
TASIUP 





Variable Description 


Initial flow (TCFS). 

Limit code for LIM (1-I, 2-d, 3-L, 
4-M ,5-P ,6-P*,7-None). 

Flow for second quarter-month (actual) 
(optional) (TCFS). 

Limit code for LIM1 (see IE). 

Sum of dissipitation (TCFS). 

Amount which running 12-month sum of 
Lake Superior plus Lakes 
Michigan-Huron net basin supplies 
must be less than in order to 
decrease the Lake Ontario outflow 
(HCFS). 

Amount which running 12-month sum of 
Lake Superior plus Lakes 
Michigan-Huron net basin supplies 
must be greater than in order to 
increase the Lake Ontario outflow 
(HCFS). 

Amount which running sum of Lake Erie 
48 quarter-month net basin supplies 
must be less than in order to 
decrease the Lake Ontario outflow 
(HCFS). 

Amount which running sum of Lake Erie 
48 quarter-month net basin supplies 
must be greater than in order to 
increase the Lake Ontario outflows 
(HCFS). 

Amount which running sum of Lake 
Ontario 48 quarter-month net basin 
supplies must be less than in order 
to decrease the Lake Ontario 
outflow (HCFS). 

Amount which running sum of Lake 
Ontario 48 quarter-month net basin 
supplies must be greater than in 
order to increase the Lake Ontario 

outflow (HCFS). 

Amount by which Lake Ontario outflow 
is to be increased (TCFS). 

Number which adjusted supply indicator 
must be greater than to increase 

the Lake Ontario outflow. 


Variable 
Columns Format Name Variable Description 


1-60 1215 ION Amount by which Lake Ontario outflows 
are to be decreased (TCFS). 

61-65 15 IASIDN Amount which adjusted supply indicator 
must be less than in order to flow 
less water from Lake Ontario 
(HCFS). 

Net basin supplies for Lake Superior 
for year previous to initial year 
(HCFS). 

Net basin supplies for Lakes 
Michigan-Huron for year previous to 
initial year (HCFS). 

Net basin supplies for Lake Erie for 
year previous to initial year 
HC). 

Net basin supplies for Lake Ontario 
for year previous to initial year 
(HCFS). 


2.3 Output Data 

The autput consists of printed tabular displays of projected stages and 
flows on each of the lakes including pertinent associated data items. In 
addition there are 19 tape outputs. 


2.3.1 Printed Output 


The printed output consists of: 1) the Lake Superior printout; 2) the 
middle lakes printout; and 3) the Lake Ontario printout. 


The Lake Superior printout contains four years of record per page. The 
record contains the monthly data for: 
The Lake Superior beginning-of-month stages (in feet-IGLD); 
The 1949 Rule flow (in TCFS); 
The Plan 1977 flow (in TCFS); 
The Compensating Works gate setting; 
tt ose flow (gate setting converted to Lake Superior outflow) 
in TCFS); 
The net total supplies for Lake Superior (in TCFS); 
The storage on Lake Superior (in TCFS and feet); 
The monthly mean stages on Lake Superior (in feet-IGLD); 
The flow retardation values for the St. Clair, Detroit and Niagara 
Rivers (in HCFS); and 
The 12-month running sums of the Lakes Superior, Michigan-Huron and 
Erie net basin supplies. 






























The middle lakes printout contains a half-year of record per page. The 
record consists of quarter-monthly data for each of Lakes Michigan-Huron, St. 
Clair and Erie. This quarter-monthly data includes the following (for each of 
the three lakes). 

1. The net total supplies for each lake (in TCFS). 

2. The quarter-monthly mean stage for each lake (in feet-IGLD). 
j 





3. The quarter-monthly mean outflow from each lake (in TCFS). 
4. The end-of-quarter month stage for each lake (in feet-IGLD). 





The Lake Ontario printout contains a half-year of record per page on a 
quarter-monthly basis. This record consists of the following. 
| 1. The maximum and minimum Plan 1958-D flows (in TCFS). 
| 
{ 
| 








2. The Lake Ontario net basin supply (in TCFS). 
3. 16.5 X the weighted supply indicator (in TCFS). 
4. The weighted supply indicator (in TCFS). 
5. The normal supply indicator (in TCFS). 
6. The adjusted supply indicator (in TCFS). 
' 7. The difference in supply indicators (in TCFS). 
; 8. The basic rule curve outflow (in TCFS). 
9. The seasonally adjusted outflow (in TCFS). 
10. The applicable limitations (in TCFS). 
11. The computed Lake Ontario outflow (in TCFS). 
12. Change in storage on Lake Ontario (in TCFS and feet). 
13. The end-of-quarter-month levels on Lake Ontario, for Plan 1958-D flow 
and for actual flow. 





2.3.2 Tape Output 
The 19 files output on tape are listed below. 


TAPE CONTENT 
80 U.S. Slip gauge monthly mean stages 

(feet-IGLD). 

Lake Superior monthly mean stages 
(feet-IGLD). 

Lake Superior monthly mean outflows (TCFS). 

Lakes Michigan-Huron monthly mean stages 
(feet-IGLD). 

Lakes Michigan-Huron monthly mean out fl ows 
(TCFS). 

Lake Erie monthly mean stages (feet~-IGLD). 

Lake Erie monthly mean outflows (TCFS). 

Lake Ontario monthly mean stages (feet-IGLD). 

Lake Ontario monthly mean outflows (TCFS). 

Lake St. Louis quarter-monthly mean stages 
(feet-IGLD). 

Lake St. Clair beginning-of-month stages 
(feet-IGLD). 

Lake St. Clair monthly mean stages 
(feet-IGD). 

Lake St. Clair monthly mean outflows (TCFS). 

Lake Superior beginning-of-month stages 
(feet-IGLD). 

Lakes Michigan-Huron beginning-of-month 
stages (feet-IGLD). 

Lake Erie beginning-of-month stages 
(feet-IGLD). 

Lake Ontario end-of-month stages (feet-IGLC). 

Lake Superior monthly mean stages 
(feet-1GLD). 

Lakes Michigan-Huron monthly mean stages 
(feet-IGD). 

Lake Erie monthly mean stages (feet-IGLD). 

Lake Ontario monthly mean stages (feet-I@0). 

Lake Superior monthly mean outflows (TCFS). 

Lakes Michigan-Huron monthly mean out fl ows 
(TCFS). 

Lake Erie monthly mean outflows (TCFS). 

Lake Ontario monthly mean outflows (TCFS). 

Lake Erie quarter-month outflows (TCFS). 

Lake Ontario quarter-month outflows (TCFS). 

Lake St. Louis quarter-month outflows (TCFS). 

Lake Ontario beginning-of-quarter-month 
stages (feet-IGD). 

Lake Ontario quarter-monthly mean stages 
(feet-IGLD). 

Lake Superior Compensating Works gate 
settings. 
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3.1 List 


BOME - 
BOMEJ - 


BOMH - 
BOMHJ - 
BOMO - 
BOMS - 
BOMSC - 
BOMSCJ - 
BOMSJ - 


Section 3 
LIST OF VARIABLES 


of Variables -- Main Program 
-A- 


Input - array containing additions to Lake Erie outflows (TCFS). 

Input - weighting factor in Lake Superior Basic Rule flaw (TCFS). 

Array of Lake Erie monthly mean levels (feet-IGLD). 

Array of Lakes Michigan-Huron monthly mean levels (feet-IGLD). 

Array of Lake St. Clair monthly mean levels (feet-IGLD). 

First mltiplicitive constant in gate-discharge equation for Lake 
Superior outflows. 

Second multiplicitive constant in gate-discharge equation for Lake 
Superior outflows. 

Additive constant in gate-discharge equation for Lake Superior 
out fl ows. 


Array containing Lake Erie beginning-of-month levels (feet-IGLD). 
Lake Erie beginning-of-month level used in forecast to determine Lake 
Superior gate setting (feet-IGLD). 


Array containing Lakes Michigan-Huron beginning-of-month levels 
(feet-IGLD). 

Lakes Michigan-Huron beginning-of-month level used in forecast to 
determine the Lake Superior gate setting (feet-IGLD). 

Array containing the Lake Ontario beginning-of-month levels 

Array containing the Lake Superior beginning-of-month levels 
(feet-IGLD). 

Array containing the Lake St. Clair beginning-of-month levels 
(feet-IGLD). 

Lake St. Clair beginning-of-month level used in forecast to determine 
gate setting (feet-IGD). 

Lake Superior beginning-of-month level used in forecast to determine 
Lake Superior gate setting (feet-IGLD). 


ee 


aE iid containing the additional Chicago Diversion flows 

HCFS). 

Adjusted Chicago Diversion outflow used in regulation (HCFS). 

a ee Chicago Diversion outflow to which adjustments are made 
HCFS). 

Input - Lake St. Clair end-of-period elevation (feet-IGLD). 

Input = allowable change in Lake Supertor flows per month (TCFS). 
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EEOP - 
EICE - 


FFOUR - 
FNIA - 
FTHREE - 


F3 - 
F4 - 


HEOP - 
HK - 


IBK - 
ICE - 


ICEC 
ICED 


ICEE 
ICEM 
ICEN 


ICESC - 


a 


Lake Erie end-of-period elevation (feet-IGLD). 
Array of flow retardation values for Niagara River (HCFS). 


ie fs 


Lake Superior outflow from basic equation after outflow 
limitations; used in forecast (TCFS). 

Sum total flow to be added (ADD) or subtracted (SUB) to/from Lake 
Erie outflow. 

a had outflow calculated by basic equation; used in forecast 
TCFS). 

Lake Superior outflow calculated by basic equation (TCFS). 

Lake Superior outflow from basic equation after outflow limitations 
are applied {TCFS). 


a Hid 


Zero, used in subroutine CANADA (can be changed to add or subtract 
to/from Lakes Michigan-Huron Elevation for changing conditions, 
i.e., 1983 to 1968 etc.). 

Lakes Michigan-Huron end-of-period elevation (Feet-IGD). 

Beginning-of-month Lakes Michigan-Huron stage, minus Lakes 
Michigan-Huron mean stage divided by Lakes Michigan-Huron standard 
deviation; used in both actual and forecast computations. 


Sfis 


Monthly index for computation of U.S. Slip gauge elevation by 
successive approximation. 

Blank character, used for blank characters put out at end-of-file. 

Array of ice retardation values for Lakes Superior, Michigan-Huron, 
St. Clair and Erie (HCFS). 

Array of ice retardation values for Lake St. Clair for period of 
record (HCFS). 

Input parameter: 1-for ice retardation in Detroit River (Lake St. 
Clair outflow); O-for no retardation. The retardation value is 
multiplied by this value and increased by a factor of N if ICED = 
N or reduced by a factor of Nif ICED =-N. 

Ice retardation values for Lake Erie for printout (fixed point to 
facilitate printing). 

Array of ice retardation values for Lakes Michigan-Huron for period 
of record. 

Input parameter: l-for ice retardation Niagara River (Lake Erie 
outflow); O-for no retardation. Since the ice retardation values 
are mltiplied by this integer, it is possible to increase the 
retardation value and/or change its sign). 

Input parameter: 1-for ice retardation in St. Clair River (Lakes 
Michigan-Huron outflow); O-for no ice retardation. Since the 
retardation values are multiplied by this integer, it is possible 
to increase the retardation value and/or change its sign). 
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IEFLW - 
IEY - 
IFO - 


IFSG - 
IFWSG - 


IGATE - 


IIG - 
III - 


IJ - 
IJK - 


INS = 


INSE - 
INS2 - 


10S - 
ISYM - 


Array of Lake Erie quarter-month flows used as inflow to Lake 
Ontario (HCFS). These values are saved for running alternate 
scenarios on Lake Ontario utilizing the same upper lakes data). 

Input - ending year. 

Input - deviation code for Lake Ontario routing Plan 1958-D: 0-for 
no deviation; 1l-for deviation. 

Input - gate setting indication. Use 10411 for Plan 1977. 

Internal variable used to determine starting and ending summer gate- 
setting months. 

Array containing gate setting for Lake Superior with 0 = no gates, 
1 = half gate, and N = N-l gates for N greater than 1. 

Output variable: either 1/2 or N/A, as required. 

Subroutine argument; used in call of subroutine PTITLE; 1l-returned 
for no end-of-file encountered. 

Do-Loop index; used to read last fifteen 'cards' of basic data; read 
in for forecast. 

Yearly mean stage index; used to increase or decrease Chicago 
Diversion: l-yearly mean stage larger than target stage (applies 
Chicago Diversion correction); O-yearly mean stage not larger than 
target stage. No Chicago Diversion correction applied). 

Same as IJK, except used in gate-setting forecast. 

Input variable: 1-for Lakes Michigan-Huron net basin supplies 
utilized in running sum (called NBH), added to Lake Erie outflow; 
0-for no utilization. 

Input variable: 1l- for Lakes Michigan-Huron net basin supplies 
utilized in running sum (called NBHE), subtracted from Lake Erie 
outflow; 0-for no utilization. 

Input parameter: 1l-for Lakes Michigan-Huron net basin supplies 
being used in NBH2 (additional Chicago Diversion water index); 
ret Lakes Michigan-Huron net basin supplies not being used in 
NBH2. 

Input parameter: 1-for Lake Superior net basin supplies being used 
in NBH (additional Lake Erie flow indication); 0-for Lake 
Superior net basin supplies not being used in NBH. 

Input parameter: l- For Lake Superior net basin supplies. 

Input parameter: 1-For Lake Superior net basin supplies being used 
in NBH2 (additional Chicago Diversion water index); 0-for Lake 
Superior net basin supplies not being used in NBH2. 

Output - Lake Superior outflow (TCFS). 

Array containing the symbols N/A and 1/2. 


ITITLE - Array containing the title of the run. 


IY - 


IYA - 
IYC - 
IYEAR - 
IYR - 


IZ - 


Number of years from beginning year necessary to read supplies for 
forecast, in main body of program for current year. 

Number of years processed. Used for printout of Lake Supertor. 

Index-year, do-loop index. 

Year. Read in from basic data on Tape 66. 

suet of years processed plus 1. Used in forecast for lower lakes 
supplies. 

Index, running from begining-year to ending-year, to read in 
forecasted supplies. 
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MAXNCD - 
MM - 


oe 


Index, month indicator. 

Lake Superior outflow indicator: 1-for December, 2-for 
January-April, and 3-otherwise (input parameter for QSUP). 

Index, month indicator for forecast. 

Index, quarter-month indicator for forecast and month loop index 
for summing stages of forecast. 


a a 


First month of accounting year for controlling Chicago Diversion by 
stage. 


2h & 


Index, used to calculate Lake Superior outflow by successive 
approximation loop (four times). 

Long Lac-Ogoki Diversion reduction amount (HCFS). 

Index: 1-for Lake Superior net basin supplies being used in NBH3 
(Long Lac-Ogoki Diversion on/off parameter); 0-for Lake Superior 
net basin supplies not being used in NBH3. 

Index: l-for Lakes Michigan-Huron net basin supplies being used in 
NBH3; ands Lakes Michigan-Huron net basin supplies not being used 
in NBH3. 


Maximum number of title cards allowed (currently 10). 
Starting month of summer forecast period. 


Array containing main program title. 

Monthly average of 12-month running sum of Lakes Michigan-Huron 
supplies (HCFS). 

Input - number which running 12-month sum of Lake Erie and/or Lakes 
Michigan-Huron net basin supplies must be less than in order to 
reduce the Lake Erie outflow. 

ners). average Lake Erie net basin supplies for current month 
HCFS). 

Input - number which running 12-month sum of Lake Erie and/or Lakes 
Michigan-Huron net basin supplies must be greater than or equal 
to in order to increase the Lake Erie outflow. 

Input - number which running 12-month sum of Lake Superior and/or 
Lakes Michigan-Huron net basin supplies must be less than in order 
to increase the Lake Erie outflow. 

Input - number which running 12-month sum of Lake Superior and/or 
Lakes Michigan-Huron net basin supplies must be greater than or 
equal to in order to increase the Lake Erie outflow. 

Input - number which running 12-month sum of Lake Superior and/or 
Lakes Michigan-Huron net basin supplies must be Tess than in order 
to reduce the Lake Erie outflow. 





NBH2 - Input - number which running 12-month sum of Lake Superior and/or 
Lakes Michigan-Huron net basin supplies must equal or exceed to 
flow additional water in the Chicago Diversion. 

NBH3 - Input - running 12-month sum of Lake Superior plus Lakes 
Michigan-Huron net basin supplies which must be equalled or 
exceeded to flow less water from Long Lac-Ogoki Diversion. 

NBMH - Input array containing Lakes Michigan-Huron net basin supplies for 
year previous to beginning year (HCFS). 

NBR - Input array containing Lake Erie net basin supplies for year 
previous to beginning year (HCFS). 

NBS - Array containing net basin supplies for Lakes Superior, 
Michigan-Huron, and St. Clair (HCFS). 

NBSE - Array containing Lake Erie quarter-month net basin supplies for each 
year used in forecast. 

NBSEQ - Array - common - containing Lake Erie quarter-month net basin 
supplies for current year. 

NBSJ - Lake Superior net basin supplies used in forecast (HCFS). 

NBSM - Array containing Lakes Michigan-Huron net basin supplies (HCFS). 

NBSS - Array containing Lake Superior net basin supplies. (HCFS). 

NBSSC - Array containing Lake St. Clair net basin supplies (HCFS). 

NBSS2 - Input array containing net basin supplies for Lake Superior for 
year previous to beginning year. 

NG - Total number of months in forecast period. 

NGATE - Output - nuber of gates used to determine Lake Superior outflow. 

NG1 - Total mumper or gates in forecast period with zero gates = 0, 1/2 
gate = 1 and N ga*2s = N+ 1 for N between 1 and 16 gates, 
inclusive. 

NHC - Parameter used in PTITLE - output - number of header cards. 

NH5 - Array containing 5% Lakes Michigan-Huron net basin supplies for 
forecast (HCFS). 

NH50 - Array containing 50% Lakes Michigan-Huron net basin supplies for 
forecast (HCFS). 

NN - Ending month of forecast period. 

NS5 - Array containing 5% Lake Superior net basin supplies used in 
forecast (HCFS). 

NS50 - Array containing 50% Lake Superior net basin supplies used in 
forecast (HCFS). 

N12E - Running 12-month sum of Lake Erie net basin supplies (HCFS). 

N12EJ - ge hs eae sum of Lake Erie forecast net basin supplies 

HCFS). 
N12MH - gat as sum of Lakes Michigan-Huron net basin suppl ies 
HCFS). 

N12MHJ - Running 12-month sum of Lakes Michigan-Huron forecast net basin 
supplies (HCFS) 

N12S - Running 12-month sum of Lake Superfor net basin supplies (HCFS). 

N12SJ - sa he aad sum of Lake Superior forecast net basin supplies 

CFS). 
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cite 


Array - common - Lake Erie monthly mean outflows (TCFS). 

Long Lac Ogoki Diversion outflow in forecast (HCFS). 

Array - common - Lakes Michigan-Huron monthly mean outflows (TCFS). 

Array - common - Lake Ontario monthly mean outflow (TCFS). 

Array containing Lake Superior outflows (part of argument of 
CANADA). 

Array containing Lake St. Clair monthly mean outflows (TCFS). 

ge is outflow for forecast month (part of argument of 
CANADA). 

Lake Superior outflow from gate setting (also parameter of QSUP). 


me ee 


Array containing Lakes Michigan-Huron target stages (feet-IGLD). 

Array containing Lakes Michigan-Huron standard deviations (feet). 

Array containing monthly mean Lakes Michigan-Huron levels for year 
previous to current year (feet-IGD). 

Array containing the Lakes Michigan-Huron monthly mean levels for 
year previous to current year; used in forecast (feet-IGD). 

Array containing Lake Superior target stages (feet-IGLD). 

Array containing Lake Superior standard deviations (feet). 


-Q- 


Internal parameter for computing U.S. Slip gauge water level (Lake 
Superior outflow divided by 1.283). 

Input array containing Lake Superior target flows (TCFS). 

Parameter in call to CANADA for monthly mean level for Lakes 
Michigan-Huron, used in forecast. 


eee 


Internal parameter, currently 0.0, used to change the level of Lake 
St. Clair by a constant for running different conditions; i.e. 
1924 or 1978. 

Array containing Lake Superior beginning-of-period stage 
(feet-IGLD). 

Array containing Lake Superior ice-retardation parameters (HCFS). 

Intermediate parameter used in calculating Lake Superior out flaw. 
Equal to adjusted normal level of Lakes Michigan-Huron. 

Lake Superior net total supply (Lake Superior net basin supply plus 
Long Lac-Ogoki Diversion). 

Array containing Lake Superior outflows (TCFS-rounded to nearest 10 


C e 

Array containing Lake Superior outflows (TCFS-rounded to nearest 10 
CFS), used in forecast. 

Array - common - Lake Erie monthly mean levels (feet-IGLD). 

Array - common - Lakes Michigan-Huron monthly mean levels 

Array - common ~ Lake Ontario monthly mean levels (feet-IGLD). 

Output - difference between Lake Superior net total supply and Lake 
Superior outflow (TCFS). 
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STS - 
STSC - 
STSJ - 
SUB - 
SUM - 


SUPLIM - 


TEMP - 
TEMPL - 


TEMPS - 
TEMPSO - 


TEMPU - 
TMMS - 


Array - common - Lake Superior monthly mean levels (feet-IGLD). 

Array - common - Lake St. Clair monthly mean levels (feet-IGLD). 

Lake Superior end-of-month level used in forecast (feet-IGLD). 

Array containing subtactions to Lake Erie outflows (TCFS). 

Average Lakes Michigan-Huron monthly mean levels for accounting year 
(feet-IGLD). 

Array containing Lake Superior lower flow limit (TCFS). 


es ae 


Input array containing target stages for increasing Chicago 
Diversion outflows. When Lakes Michigan-Huron beginning-of-month 
level is greater than the corresponding target stage, the Chicago 
Diversion outflow may be increased. 

Intermediate parameter. 

Change limit in Lake Superior lower flow limit (previous month flow 
minus allowable change (CHNG)) (TCFS). 

Lake Superior previous month outflow, used in forecast (TCFS). 

Input - Lake Superior monthly mean flow for the month previous to 
the starting month (TCFS). 

Upper Lake Superior outflow change limit (previous month flow plus 
allowable change (CHNG) (TCFS). 

Input - target level for increasing Chicago Diversion for specific 
eer if previous year average leve) exceeds or is equal to this 

evel. 


= 
Array containing the U.S. Slip gauge levels (feet-IGLD). 


Initial guess of U.S. Slip stage (feet-IGLD). 


Pe ae 
Additional Chicago Diversion index: l-when the running sums or 
supplies exceed targets; O-otherwise. Also used as intermediate 
parameter in calculating the U.S. Slip stage. 
eos 
Additional Chicago Diversion index: l-when the Lakes Michigan-Huron 
beginning-of-month level exceeds monthly target level; 
O-otherwise. Also used as intermediate parameter in calculating 
the U.S. Slip stage. 
ay se 


0.0; used to write final record for middle lakes print file. 





3.2 List of Variables -- Subroutines 


The main program P4827 uses 11 subroutines. These are listed and briefly 
described below. 


3.2.1 Subroutine ABSICE 


Purpose. To compute absolute ice where ice is a floating point array. 
Using intrinsic function ABS, returns absolute value of ice. 


Argument list. ABSICE (AICE). 


Description of variables. 


AICE - Input and output - vector, length 12, of which the absolute value is 
taken. 
I - Loop index. 


3.2.2 Subroutine QSUP 
Purpose. Compute Lake Superior gate settings and flow. 
Arguments. QSUP(JJ, REQFL, SBEG, OUT, X1, Y1, Z1, N, M,). 


Description of variables. 


DIF - Internal - difference between regulated flow and gate flow. 

JJ - Input - month index: 1=Dec; 2=Jan to Apr; 3=May to Nov. 

M - Input - index for calculating gate setting; 0 for returning constants 
for given gate settings, 1 for calculating gate setting. 

N- Input and output - number of gates; 1=1/2 gate, 2=1 gate.... 17=16 
gates. 

OUT - Qutput - Lake Superior outflow for gate setting (TCFS). 

REQFL - Input - Lake Superior regulated flow (TCFS). 

SBEG - Input - Lake Superior beginning-of-month level (feet). 

X- Array containing conversion constant for gate setting. 

Xl - Output - feet to TCFS conversion constant for gate setting. 

Y- Array containing stage factor constants for gate equations. 

Yl - Output - stage factor constant for gate setting. 

Z- Array containing mean bottom constants for gate equations. 

Zl - Output - mean bottam constant for gate setting. 


3.2.3 Subroutine CANADP 


Purpose. Subroutine to route middle lakes, i.e. campute monthly stages 
and out Tow for Lakes Michigan-Huron, St. Clair and Erie, using 1962-1968 
conditions. 








me eee 
ES AR nr 


me 


Arguments. 


ADD - Common - 
AVEE = Common - 
AVEMH - Common - 
AVESC - Common - 
BOME - Common - 
BOHM - Common - 
BOMO - Common - 
BOMS - Common - 
C - Internal 


D- Internal Parameter - Lake St. Clair stage - 543.81 feet used in Lake 
St. Clair flow equation. . 


EOPE - Input and sutput - end-of-period level Lake Erie (feet-IGLD). 
EQPMH - Input and output - end-of-period level Lakes Michigan-Huron (feet- 


IGLD). 


EOPSC - Input and output - end-of-period level Lake St. Clair (feet). 


CANADP (IYR, K, EOPMH, EOPSC, EOPE, SMH, RISC, SLSC, RID, RIN 
XAV3, XAV6, XAV9) 


Description of variables. 


array containing additions to Lake Erie outflow. 
array of Lake Erie monthly mean levels. 

array of Lakes Michigan-Huron monthly mean levels. 
array of Lake St. Clair monthly mean levels. 


array of Lake Erie beginning-of-month levels. 

array of Lakes Michigan-Huron beginning-of-month levels. 
array of Lake Ontario beginning-of-month levels. 

array of Lake Superior beginning-of-month levels. 


2 
Lakes Michigan-Huron outflow equation. 


CD - Common - array of additional Chicago Diversion flows (TCFS). 
COP - Common - basic Chicago Diversion flow (TCFS). 


EQNBS - Internal - 


HE - Internal 
HMH - = Internal 
HMH1 - Internal 


(used in 


HSC - = _ Internal 
HSCl1 - Internal 


flow equation). 


I - Index of loop to calculate average Lakes Michigan-Huron stage for 
accounting year. 

IBY - Common - Beginning year of regulation. 

IES - Common - Lake Erie quarter month outflow for running other Lake 
Ontario plans (HCFS). 

IEY - Common - ending year of regulation. 


IFOP - Common - Deviation code for Lake Ontario with 1958-D, l=yes 0=no. 

IJK - Common - yearly mean stage index; used to increase or decrease Chicago 
Diversion: l-yearly mean stage larger than target stage (applies 
Chicago Diversion correction); O-yearly mean stage not larger than 
target stage (no Chicago Diversion correction applied). 

INM - Common - 1-for Lakes Michigan-Huron NBS utilized in running sum (NBH) 
(added to Lake Erie outflow); O-for no utilization. 

INME - Common - l-for Lakes Michigan-Huron supplies in NBHE (subtract from 
Lake Erie outflow); O0-otherwise. 






Lake Erie quarter month NBS (TCFS). 


average Lake Erie level for period. 

average Lakes Michigan-Huron level for period. 

average Lakes Michigan-Huron level for period, adjusted 
flow equation). 

average Lake St. Clair level for period. 

average Lake St. Clair level for period, adjusted (used in 
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INM2 - Common - l-for Lakes Michigan-Huron supplies in NBH2 (add to Chicago 
Diversion); O-otherwise. 

INS - Common - 1l-for Lake Superior supplies in NBH (add to Chicago 
Diversion); 0-otherwise. 

INSE = Common - l-for Lake Superior supplies in NBHE (decreases Lake Erie 
outflow); O-otherwise. 

INS2 - Common - 1-for Lake Superior supplies in NBH2 (additional Chicago 
Diversion); O-otherwise. 

IY - Common - current regulation year. 

IYR - Input - current regulation year. 


J- Index of tenth-of-quarter loop. 
JDX - Internal - end-of-month index: 1-for end-of-month; 2-otherwise. 


K - Input - current regulation month. 
Kl - Common - first month of accounting year for controlling Chicago 
Diversion by stage. 


L - Index of successive approximation loop. 
LA - Quarter-of-year index, used for converting Lake Erie supplies. 
LL - Index for loop converting Lake Erie supplies. 


Ml - Internal - present quarter-of-month + 1. 
M2 - Internal - present quarter-of-year. 


NBEE ~ Common - number which running twelve-month sum of Lakes Erie and/or 
Michigan-Huron NBS must be < to reduce Lake Erie out fl mw. 
NBE1 ~ Common = erage Lake Erie NBS for current month. 


NBE2 - Common - number which running twelve-month sum of Lakes Erie and/or 
Michigan-Huron must be > to increase Lake Erie outflow. 

NBE3 = Common - number which running twelve-month sum of Lake Erie NBS must 
be > to increase Chicago Diversion. 

NBH - Common - number which running twelve-month sum of Lakes Superior 
and/or Michigan-Huron NBS must be > to increase Lake Erie outflow. 

NBHE - Common - number which running twelve-month sum of Lakes Superior 
and/or Michigan-Huron NBS must be < to reduce Lake Erie outflow. 

NBH2 = Common - number which running twelve-month sum of Lakes Superior 
and/or Michigan-Huron NBS must > to increase Lake Erie flow. 

NBS - Common - array of Lakes Superior, Michigan-Huron and St. Clair NBS. 

NBSEQ - Common - array of Lake Erie quarter-month NBS for current regulation 
year. 

N12E - Common - running twelve-month sum of Lake Erie NBS (HCFS). 

N12MH - Common - running twelve-month sum of Lakes Michigan-Huron NBS (HCFS). 

N12S - Common - running twelve-month sum of Lake Superior NBS (HCFS). 


OE - Internal - outflow from Lake Erie. 

OM - Common - array of outflows for Lakes Michigan Huron, monthly means for 
current regulation year. 

OMH - Internal - outflow Lakes Michigan-Huron. 

00 - Common - array of Lake Ontario monthly mean outflow for current 
regulation year. 

OSC - Internal - outflow of Lake St. Clair. 
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Common - array of Lakes Michigan-Huron monthly mean levels for 12 
months previous to current regulation month. 

Tape 60 input - Lakes Michigan-Huron quarter-month NTS (TCFS). 

Tape 60 input - Lakes Michigan-Huron quarter-month stage (feet). 

Tape 60 input - Lakes Michigan-Huron quarter-month flow (TCFS). 

Tape 60 input - Lake St. Clair quarter-month NTS (TCFS). 

Tape 60 input - Lake St. Clair quarter-month stage (feet). 

Tape 60 input - Lake St. Clair quarter-month flow (TCFS). 

Tape 60 input - Lake Erie quarter-month NTS (TCFS). 

Tape 60 input - Lake Erie quarter-month stage (feet). 

Tape 60 input - Lake Erie quarter-month flow (TCFS). 

Internal - quarter-month sum Lakes Michigan-Huron NBS (TCFS). 

Internal - quarter-month sum Lakes Michigan-Huron stage (feet). 

Internal - quarter-month sum Lakes Michigan-Huron flow (TCFS). 

Internal - quarter-month sum Lake St. Clair NBS (TCFS). 

Internal - quarter-month sum Lake St. Clair stage (feet). 

Internal - quarter-month sum Lake St. Clair flow (TCFS). 


Internal - quarter-month sum Lake Erie NBS (TCFS). 
Internal - quarter-month sum Lake Erie stage (feet). 
Internal - quarter-month sum Lake Erie flow (TCFS). 


Input - retardation of Lake St. Clair outflow. 
Input - retardation of Lake Erie outflow. 
Input - retardation of Lakes Michigan-Huron outflow. 


Common - array of Lake Superior monthly mean levels, current 
regulation year. 
Input - Lake St. Clair NBS (TCFS). 
Input - Lakes Michigan-Huron NBS (TCFS). 
Internal - Lake Erie quarter-month NTS. 
SNTMH - Internal - Lakes Michigan-Huron quarter-month NTS minus Chicago 
Diverion. 
SNTSC - Internal - Lake St. Clair quarter-month NTS. 
SO - Common - Lake Superior flows, current regulation year. 
SOT - Common - Lake Ontario levels, current regulation year. 
STE - Common - Lake Erie levels, current regulation year. 
STMH - Common - Lakes Michigan-Huron levels, current regulation year. 
SUB - Common - amount deducted from Lake Erie outflow. 
SUM - Internal - average of Lakes Michtgan-Huron monthly mean levels for 
preyious accounting year. 


TBOMS - Common - target levels for Chicago Diversion. Lakes Michigan-Huron 
beginning-of-month levels must exceed coordinated level to increase 
diversion. 

TMMS - Common - target level for Chicago Diversion. Chicago Diversion 
increases if yearly average Lakes Michigan-Huron monthly mean for 
level accounting year exceeds this. 


Internal - Chicago Diversion supply index. When twelve-month 
moving sum supplies exceed targets x=l (Chicago Diversion 
increased); x=0 otherwise (no increase). 

Be Tape 60 input - monthly mean Lakes Michigan-Huron NTS 
TCFS). 
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XAV2 - Internal, Tape 60, input - monthly mean Lakes Michigan-Huron stage 
(feet-IGLD). 

XAV3 - eS). Tape 60, input - monthly mean Lakes Michigan-Huron fl ow 
TCFS). 

XAV4 - Internal, Tape 60, input - monthly mean Lake St. Clair NTS (TCFS). 

XAV5 - Internal, Tape 60, input - monthly mean Lake St. Clair stage 
(feet-IGLD). 


XAV6 - Internal, Tape 60, input - monthly mean Lake St. Clair flow (TCFS). 
XAV7 - Internal, Tape 60, input - monthly mean Lake Erie NTS (TCFS). 

XAV8 - Internal, Tape 60, input - monthly mean Lake Erie stage (feet-IGLD). 
XAV9 - Internal, Tape 60, input - monthly mean Lake Erie flow (TCFS). 

XHE - Internal - Lake Erie beginning-of-period levels (feet-IGLD). 

XHMH - Internal - Lakes Michigan-Huron beginning-of-period levels 

XHSC Internal - Lake St. Clair beginning-of-period levels (feet-IGLD). 


XSUML ~ Internal - ra Michigan-Huron monthly sum of tenth-of-quarter-month 
NTS (TCFS). 

XSUM2 - Internal - Lakes Michigan-Huron monthly sum of tenth-of-quarter-month 
stages (feet-IGLD). 

XSUM3 ~ Internal - Lakes Michigan-Huron monthly sum of tenth-of-quarter-month 
flows (TCFS). 

XSUM4 - Internal - Lake St. Clair monthly sum of tenth-of-quarter-month NTS 


(TCFS). 

XSUM5 ~ Internal - Lake St. Clair monthly sum of tenth-of-quarter-month stages 
(feet-IGLD). 

XSUM6 - Beat - Lake St. Clair monthly sum of tenth-of-quarter-month flows 
TCFS). 


XSUM7 - Internal - Lake Erie monthly sum of tenth-of-quarter-month NTS (TCFS). 
XSUM8 - Internal - Lake Erie monthly sum of tenth-of-quarter-month stages 


(feet-IGLD). 

XSUM9 - Internal - Lake Erie monthly sum of tenth-of-quarter-month fl ows 
(TCFS). 

XX - Common - Lake Erie monthly mean flow, for current regulation year. 

Y- Internal - 1-for Lakes Michigan~Huron beginning-of-month level 
exceeding Chicago Diversion increase target level for month, 
0-otherwise. 


3.2.4 Subroutine CANADA 


Purpose. Subroutine used to canpute the average gate setting for the Lake 
Superior Compensating Works. 


Arguments. CANADA (A, B, IYR, K, K1, EOPMH, EOPSC, EOPE, SMH, RISC, SLSC, 
ID, SE, RIN, NBSE, CD, FNIA, QSUM). 


Description of variables. 


A - Input - number added to Lakes Michigan-Huron stage to change 
conditions. 


B- Input - number added to Lake St. Clair stage to change conditions. 


C - Internal parameter - used in Lakes Michigan-Huron outflow equation: 
[(M-H stage + St. C stage)/2]-543.40. 
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cD- Input - Chicago Diversion. 


D- Internal parameter - used in Lake St. Clair outflow equation: 
St. C stage - 543.81. 


EOPE - Input and output - end-of-period Lake Erie level. 
| EOPMH - Input and output - end-of-period Lakes Michigan-Huron level. 
EDPSC - Input and output - end-of-period Lake St. Clair level. 
| 


FNIA - Input - additional Lake Erie outflow. 


HE - Internal - average Lake Erie level for period. 

HMH - = Internal - average Lakes Michigan-Huron stage for peiord. 

HMHl - Internal - average Lakes Michigan-Huron stage for period, corrected 
(used in flow equation). 


- ee ee 


HSC - Internal - average Lake St. Clair elevation for period. 
HSC1 - Internal - average Lake St. Clair elevation, corrected (used in flow 
equation). 


IYR - Input - not used. 


J- Index of tenth-of-quarter loop. 
K - Input - not used. 
Kl - Input - not used. 
L - Index of successive approximation loop. 


i LL - Index of quarter-month loop. 
NBSE - Input - Lake Erie quarter-month supplies (HCFS). 


OE - Internal - outflow from Lake Erie for period (TCFS). 
OMI - Internal - outflow from Lakes Michigan-Huron for period (TCFS). 
OSC - Internal - outflow from Lake St. Clair for period (TCFS). 


QSUM - Output - Lakes Michigan-Huron monthly mean level (feet-IGLD). 
RID - Input - retardation - Lake St. Clair outflow. 


SE - Input - not used. 

SLSC - Input - Lake St. Clair NBS (TCFS). 

SMH - Input - Lakes Michigan-Huron NTS (TCFS). 

SNTE - Internal - Lake Erie NTS (TCFS). 

SNTMH - Internal - Lakes Michigan-Huron NTS minus Chicago Diversion (TCFS). 
SNTSC - Internal - Lake St. Clair NTS (TCFS). 


A NE NR * 


XHE - Internal - Lake Erie beginning-of-period level. 
XHMH - Internal - Lakes Michigan-Huron beginning-of-period level. 
XHSC - Internal - Lake St. Clair beginning-of-period level. 
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3.2.5 Subroutine R1958D 


Purpose. Compute stages and flows for Lake Ontario given beginning stage 
and net total supplies (NTS). 


Description of Parameters. 


AVEE - 
AVEMH - 
AVESC - 


BOQ - 


CON - 
CON139 - 
CONI7 - 
CON48 - 
CON6 - 
CON7 - 
CONS - 


EEl - 


EH 

£01 

E02 - 
ES] - 
HE2 

HSL 

[= 
TACUM - 
TASIDN - 
TASIUP - 


IBY - 
IC - 


ICE - 
ICESL - 
INDEX - 
INK - 
IOH - 


Common - average monthly mean stage Lake Erie (feet-IGLD). 

Common - average monthly mean stage Lakes Michigan-Huron (feet-IGLD). 
Common - average monthly mean stage Lake St. Clair (feet-IGLD). 
Internal - array with beginning-of-period levels Lake Ontario. 
Internal - 1/16.5. 
Internal - 1/139. 
Internal - 1/17. 
Internal - 1/48. 
Internal - 1/6. 
Internal - 1/7. 
Internal - 1/9. 


Common - Lake Erie beginning-of-month levels during current 
regulation year. 

Common - Lake Erie monthly mean levels during current regulation 
year. 

Common - Lakes Michigan-Huron beginning-of-month level, current 
regulation year. 

Common - Lake Ontario beginning-of-month level, current regulation 
year. 

Common - Lake Erie month mean flows, current regulation year. 

Common - Lake Superior beginning-of-month levels, current regulation 
year. 


Common - Lakes Michigan-Huron monthly mean levels, current regulation 
year. 
Lake St. Louis quarter monthly mean level, current regulation year. 


Internal - current quarter-of-month, index of loop. 

Internal - accumulative deviation from plan (TCFS) 

Input - adjusted supply indicator which must not be exceeded for Lake 
Ontario flow to be increased. 

Input - adjusted supply indicator which must be exceeded for Lake 
Ontario flow to be increased. . 

Common - beginning year of regulation. 

mL ne a input - symbol of applicable limitation (blk, I, d, 

’ > s 

Internal - Lake Ontario ice retardation, current regulation year. 

Internal - Lake St. Louis ice retardation, current regulation year. 

Input - "I" limitation index: 1-for I limit, 0-otherwise. 

Internal - change in stage (hundreths of a foot). 

Input - number which running sum of Lake Ontario 48 quarter-month NBS 
must be greater than to increase Lake Ontario flow. 
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ISA - 


ISI - 
ISIT - 
IST - 
ISTAG - 


ISTAGE - 


ISTAG] - 
ISTOR - 
ISUM - 
ITOT - 
IUP - 
IWS - 

IX - 

IXY - 

IY - 
IYEARB - 
IYEARE - 
IYR - 


JSTAGE - 
IDIS - 
IDN - 


IE - 
IEH - 


Input - number which running sum of Lake Ontario 48 quarter-month NBS 
must be less than to decrease Lake Ontario outflow. 

Input - Basic Rule Curve flow. 

Internal - adjustment in Basic Rule Curve flow. 

Input - quarter of first computations line printout. 

Input - final quarter. 

Write Tape 61 - current quarter of regulation month. 

Internal - used to determine appropriate flow limit. 

Write Tape 61 - Basic Rule Curve flow plus adjustment. 

Write Tape 61 - flow after limit applied. 

Write Tape 61 - supply indicator (TCFS). 

ai: ae 61 - seasonal adjustment to Basic Rule Curve flow 
TCFS). 

Input - previous 22 quarter supply indicators. 

Internal-initial change in supply indicators in quarters. 

Write Tape 61 - 16.5 x weighted supply indicator (TCFS). 

Internal - current beginning-of-period stage minus 24,000, used in 
basic flow computation (0.01 of a foot). 

Internal - current beginning-of-period stage minus 24,000, used in 
"L" limitation computation (0.01 of a foot). 

Internal - current end-of-period stage minus 24,000 (0.01 of a foot). 

Write Tape 61 - storage on Lake Ontario (TCFS). 

Write Tape 61 - adjustment to supply indicator (TCFS). 

Write Tape 61 - adjustment to supply indicator (TCFS). 

Input - amount Lake Ontario flow is increased by, by month (TCFS). 

Write Tape 61 - weighted supply. 

Internal - current quarter-of-year. 

Internal - current quarter-of-month. 

Common, write Tape 61 - current regulation year. 

Input - year of first line printed out. 

Input - ending year of regulation period. 

Input Tape 66 - current year. 


Write Tape 61 - end-of-quarter theoretical level (0.01 of a foot). 


Input Tape 66 - quarter-month deviation from plan, current regulation 
year. 

Input - amount/month Lake Ontario flow decreased by (TCFS). 

Input - limitation code number: l=I, 2=J, 3=L, 45M, S=P, 6=P* 7=bik. 

Input - running sum of Lake Erie 48 quarter-months NBS which must be 
exceeded for the Lake Ontario outflow to be increased. 

Input - running sum of Lake Erie 48 quarter-months NBS which must not 
be exceeded if the Lake Ontario. outflow to be decreased. 

Common - Lake Erie quarter-month flow (HCFS). 

Input - exceedance probability (0). 

Conmon - ending year of regulation. 

Input- optional limitation code for initial regulation quarter (see 
I e 

Internal - change in supply for previous nine quarters. 

Input Tape 66 - difference between Lake Ontario and Lake St. Louis 
flow for current year, by quarter. 

Common - Lake Ontario deviation code: l-deviation, O-no deviation. 
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Loop index. 

Internal - entry index. 

Internal - negative adjustment indicator: 1-for positive, 0-for 
negat ive. 

Target number that running twelve-month sum of Lakes Superior and 
SY cig NBS must be greater than to increase Lake Ontario 
out fl ow. 

Target number that running twelve-month sum of Lakes Superior and 
ee la NBS must be less than to decrease Lake Ontario 
Out fl ow. 


JSTAGE1 - Write Tape 61 - end-of-quarter actual level (0.01 of a foot). 


K - Internal - supply indicator index. 
KL - Index of loop to move supply indicators. 


L - Index of difference between supply indicators. 

LIM - Write Tape 61 - Lake Ontario theoretical flow (TCFS). 

LIM1 - Input optional - first quarter forced flow (TCFS). 

LPSTAR - Internal - used for “P® limitation (225 - Giff. in flow), 
6 


M - Internal - quarter of year in Lake St. Louis stage computation. 
MAX - Input- maxioum "P° limitation values. 
MAXE - Internal - number of appl icable maximum limit. 
MAXL - Internal - value of “L” limitation. 
MAXLIM - Internal - v’ jue of maximum limit. 
MAXP - Internal - caximm “P" limitation value for current quarter. 
MAXP1 - Internal - MAXP + previous quarter supply indicator leaximia “p™). 
MIN - Input - minim “P” limitation values. 
MINP - Internal - minimum "P” limitation value for quarter. 
MKL - Index for loop to move monthly mean levels and flows for Lake 
Ontario. 
ML - Internal - current value of "M° limitation. 
MLIM - Input - values of ” limitation by month. 
Internal - used as upper limit in quarter-month output loop. 
0B Input - beginning-of-month of first line printout. 
MOE Input - last month of regulation period. 
MON Internal - name of months. 
Internal - used to move supply indicators to new year. 
Internal - initial quarter for regulation year loop. 
M2 - Internal - final quarter for regulation year loop. 


Read from Tape 66 - Lake Ontarfo NBS for current regulation year. 
Common - Lake Erie quarter month NBS. 

Common - NBS for Lakes Superior, Michigan-Huron and St. Clair. 
Input - Lake Erie previous year NBS. 

Internal - used as lower limit in quarter-monthly output loop. 
Input - Lakes Michigan-Huron previous years NBS. 

Input - initial quarter NTS (TCFS). 

Input - Lake Ontario previous years NBS. 

Input - Lake Superior previous years NBS. 
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NSE - 
NSO - 
NSSMH - 


NTS - 
NTSO - 
NWS - 
OE2 - 
002 - 
0Qk - 
OQMS - 
OUT - 
POW - 
SE2 - 


SLOF - 
$02 - 


STAGE - 
STAG] - 
SUM - 


XNK - 


Interral - running 48 quarter-month sum of Lake Erie NBS. 

Internal - running 48 quarter-month sum of Lake Ontario NBS. 

Internal - running twelve-month sum of Lakes Superior and 
Michigan-Huron NBS. 

Internal - Lake Ontario quarter-month net total supply (TCFS). 

Internal - Lake Ontario quarter-month net total supply (HCFS). 

Write Tape 61 - normal weighted supply. 


ee ) Lake Ontario monthly mean stage for current regulation year 

TCFS). 

sue j Lake Ontario monthly mean flow for current regulation year 
TCFS). 

Interna] - Lake Ontario quarter-month flow (TCFS). 

Internal - Lake Ontario quarter-month stage (feet). 

Internal - Lake Ontario monthly mean flow for current regulation 

year (TCFS). 


Internal - power of Lake St. Louis level equation. 


igre } Lake Superior monthly mean level for current regulation year 

feet). 

Internal - Lake St. Louis quarterly flow (TCFS). 

ree ) Lake Superior monthly mean flows for current regulation year 

TCFS). 

Internal - Lake Ontario monthly mean levels for current regulation 
year (feet). 

Internal - Lake Ontario actual level (feet). 

Internal - sum and average of changes in seasonal index in 12 
quarters. 


Output - change in stage. 


XSTAGE - Output - Lake Ontario theoretical stage (feet). 
XSTAG1 - Output - Lake Ontario actual stage (feet). 


3.2.6 Subroutine PRTHCE 


Purpose. This subroutine outputs, in table form, the stages and outflows 
for Lakes 


ichigan-Huron, St. Clair, and Erie. 


Argument list. PRTHCE (NAM, MAXNCD, ITITLE, NCD). 
Description of variables. 


EOPE - 
EQPMH - 
EOPSC - 


IIT - 
ITILE - 
ITP - 
I® - 


Output - end-of-period Lake Erie elevation (feet). 
Output - end-of-period Lakes Michigan-Huron elevation (feet). 
Output - end-of-period Lake St. Clair elevation (feet). 


Parameter of CALL PTITLE (end of file indicator). 
Parameter of CALL PTITLE (title of run). 

Number of input tape = 60. 

Output - current run year. 





Index of quarter-month loop. 
Output - current run month. 

Index of half-year loop. 

Index of ninth loop for half-year. 


Parameter of CALL PTITLE (maximum number of title cards). 
Output - current run quarter. 


Parameter of CALL PTITLE (maximum number of title cards). 


Out put 
Output 


- Lakes Michigan-Huron quarter-month NTS (TCFS). 

- Lakes Michigan-Huron quarter-month mean stage (feet). 
Output - Lakes Michigan-Huron quarter-month mean flow (TCFS). 
Output - Lake St. Clair quarter-month NTS (TCFS). 

Output - Lake St. Clair quarter-month mean stage (feet). 
Qutput - Lake St. Clair quarter-month mean flow (TCFS). 

Output - Lake Erie quarter-month NTS (TCFS). 

Qutput - Lake Erie quarter-month mean stage (feet). 

Qutput - Lake Erie quarter-month mean flow (TCFS). 


Output 
Output 
Out put 
Output 


Lakes Michigan-Huron monthly NTS (TCFS). 

Lakes Michigan-Huron monthly mean stage (feet). 
Lakes Michigan-Huron monthly mean flow (TCFS). 
Lake St. Clair monthly NTS (TCFS). 


Output 
Output 
Output 
Out put 


3.2.7 Subroutine PRTONT 


Lake St. Clair monthly mean flow (TCFS). 
Lake Erie monthly mean NTS (TCFS). 

Lake Erie monthly mean stage (feet). 
Lake Erie monthly mean flow (TCFS). 


Output - Lake St. Clair monthly mean stage (feet). 


Purpose. This subroutine outputs, in tabular form, the stages and outflow 
for Lake Ontario. 


Argument list. PRTONT (NAM, MAXNCD, ITITLE, NCD). 


Description of variables. 


I - Index of month loop. 

IC - Output - letter for applicable limit (blank, I,J,L,M,P,P*) 
IFF - Output - change in supply indicator (TCFS). 

III - Parameter of CALL PTITLE (end of file index). 

IN - Input from Tape 61 - change in storage (feet) = XNK. 

IQ - Output - Basic Rule Curve flow. 

1Q2 - Output - seasonal adjusted flow (IQ + ISA) (TCFS). 

193 - Output - Lake Ontario flow after limit applied to (TCFS). 
IS - Output - supply indicator (TCFS). 

ISA - Output - seasonal adjustment to Basic Rule Curve flow (TCFS). 
IST - Output ~ 16.5 x weighted supply (TCFS). 

ISTAGE - Input from Tape 61 - theoretical stage = XSTAGE. 

ISTAG] - Input from Tape 61 - actual stage = XSTAG1. 
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ISTOR - Output - change ir storage (TCFS). 

ISUM - Output - adjustment to supply indicator (TCFS). 
ITITLE - Parameter of CALL PTITLE (title of run). 

ITOT - Output - adjusted supply indicator (TCFS). 


ITP - Number of input tape = 61. 

IWS - Output - weighted supply (TCFS). 

IY - Qutput - current run year. 

{YI - Output - current run quarter. 

LIM - Output - Lake Ontario flow after limits applied (TCFS). 
MAXNCD - Parameter of CALL PTITLE (maximum number of title cards). 
MAXP - Output - maximum "P" limit for current run quarter. 
MINP - Qutput - minimum "P" limit for current run quarter. 
MON - Output - name of run month. 

NAM Parameter of CALL PTITLE (name of program). 

NCD Parameter of CALL PTITLE (number of header cards). 


NTS - Output ~ Lake Ontario NTS (TCFS). 
NWS - Output - normal weighted supply for 1958-D. 


XNK - Output - change in stage (feet). 

XSTAGE - Output - theoretical plan stage (feet). 

XSTAG] - Output - actual plan stage (feet). 
3.2.8 Subroutine ROUND 


Purpose. This subroutine rounds floating point numbers using U.S. Lake 
Survey rounding. 


Argument list. ROUND (A, IA). 


Description of variables. 


A- Input floating point mumber. 
IA - Output fixed point number. 


3.2.9 Subroutine IABICE 


Purpose. Computes absolute value of ICE, where ICE is a fixed point 
array. 


Argument list. IABICE (ICE). 
ICE - Input, output - array containing 12 ice values. 
3.2.10 Subroutine PTITLE. 
Purpose. Subroutine to read and print date, time and titles. 
Argument list. PTITLE (111, NAM, MAXNCD, ITITLE, NCD, JuJ). 
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Description of varibles. 


[= 

IBK - 
IDAY - 
IHR - 
III - 


IMON - 
ITITLE - 
IR - 


J 
JJJ - 


M - 
MAXNCD - 


N- 
NAM - 
NCD - 


Index - print NAM loop and print ITITLE loop. 

Internal - parameter containing blank. 

Parameter in CALL STDMY (date of run). 

Parameter in CALL STDMY (hour and minute of run). 

Input - 1) read date and time; 2) read title to blank; 3) print time 
and title. 

Parameter in CALL STDMY (name of month of run). 

Output - array containing title of run. 

Parameter in CALL STDMY (year of nun). 


Index - count number of cards. 
Output - end of file index: 1-for not end of file; 2-for end of 
file. 


Start indicator of last title card. 
Input - maximum number of title cards. 


Number of character spaces title requires. 
Input - array with program name. 
Output - number of title cards. 


3.2.11 Subroutine STDMY 


Purpose. Determines date and time of run using CDC subroutines for 6000 or 


7600 machines at LBL. 


Argument list. STDMY (IDAY, JMO, IYR, JHR). 
Description of varibles. 


IDATE - 
IDAY - 
IHR - 
IMIN - 
IMON - 
ISEC - 
ITIME - 
IYR - 
JHR = 
JMO - 


Parameter of CALL DATE (contains date in form: dd/mn/yy). 
Output - day of run. 

Internal - hour of day of run. 

Internal - minutes of time beyond hour of run. 

Interna] - month of date of run. 

Internal - seconds of time beyond whole minutes of time of run. 
Parameter of CALL HOUR (time in form: hh/mn/ss). 

Output - year. 

Output - hour and minute of time of run. 

Output - name of month. 
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SAMPLE INPUT 
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——_____SAMPLE_ INPUT DATA FOR PLAN 1977 0020 RUN 
OP77 
— $00 .29600, 06599 87599 82600, 0760042600, 67600 ,82600,87600,80600,72600,53 PMS NPTY 
43 243 243 245 ,48 48 4 4a 43 43 243 4} PSSONPT7 
877, 72577 ,64577 64577, 80578, 15578 44578 63578 66578 S3578, 34578, 09577,92 PHM PPT 
4014 1.13 1e17 1.20 1.22 1,24 4,24 1,18 1,18 1,16 PHSODPTT 
5 ee be : Ps i) Of : 
5S. 55, 5S. S55. 55. SS, S55. 55, 55, 55, SUPLIMNP77 
_—-.1900 1900. 9. 9 0 RUN 
601,51 30, 200, 1006, RUN 
_____ #40, : o47, 034, eu@, =O, =o 15, = S13, 39, 26, 6, oe, CE SCE DPT? 
0 1 122055 . ME IX 
577,51. 572,50 569,95 3,2 196201968 CONDITIONS RUN 
; 1000 1000 1000 1900 1000 1900 1000 1000 1900 1on0 1000 1000 NOSS2 PML 
—___2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 2000 - NBM oe OML 
0 0 0 0 0 0 0 0 0 0 0 0 NBR 
———__1__. 01326599999 9 : 
° 0,9 0.0 0,0 0.0 0,0 0.0 0,0 0.0 0,0 0,0 0,0 0, 
——0.0.0 ..0.9000,00000,00000,00000,00000.00000, 00000, 0000, 00800, 00000, 00000, no  TAOMSNM 
000,00 o | TMMSOML 
_._578.09578,00578,00578,00578,00578, 00578, 00578, 00578, 00578,00578,00578,00 — PMMSNML 
2205599999 1 121840 2053 1 { ERPMCOML 
—____0,9 0,0 0,0 0,0 0,0 0,0 0,0 0,0 6,0 0,0 0,0 0,0 ANDNDNML 
0. 0.0 0,n 0.0 0,90 0.20 0.0 0.0 0.0 .0 0,0 0.0 SUBDMPL 
0 0 0? a8 4B 24a a8 GA 248 an 25 259 2610580 
263 265 266 268 268 268 267 26 0 0 0 0 00580 
Oo. oO =r OO CC TCC CCS AD 
0 0 0 0 0 0 0 0 0 0 0 0 0 0 
_______0 227 _ 232. 237 242 245 247 2a9 253 252 253 254 255 256 256 2560580 
256 255 252 249 247 auS Pad 244 240 239 238 asA aga 0 7) 0n5é6od 
210 207 204 in8 tan 190 393 193 193 493 
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Lake Superior Monthly Mean Elevation 

Lakes Michigan-Huron Monthly Mean Elevation 
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PLAN 25N 


Lake Superior Monthly Mean Elevation 

Lakes Michigan-Huron Monthly Mean Elevation 
Lake St. Clair Monthly Mean Elevation 

Lake Erie Monthly Mean Elevation 


Lake Ontario Monthly Mean Elevation (With Deviation) 
Category 1 


Lake Ontario Monthly Mean Elevation (With Deviation) 
Category 2 

Lake Ontario Monthly Mean Elevation (With Deviation) 
Category 3 

Lake Superior Monthly Mean Outflow 
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Category 1 
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Category 2 

Lake Ontario Monthly Mean Outflow (With Deviation) 
Category 3 

Lake St. Louis Monthly Mean Outflow (With Deviation) 
Category 1 
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Category 2 
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